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We carry out first-principles density-functional calculations of the antiferrodistortive (AFD) and 
ferroelectric (FE) soft-mode instabilities in tetragonal SrTiC>3, with the structural degrees of freedom 
treated in a classical, zero-temperature framework. In particular, we use frozen-phonon calculations 
to make a careful study of the anisotropy of the AFD and FE mode frequencies in the tetragonal 
ground state, in which an J?-point AFD soft phonon has condensed. Because of the anharmonic 
couplings, the presence of this AFD distortion substantially affects both the AFD and FE mode 
frequencies. The AFD mode is found to be softer for rotations around a perpendicular axis (E g 
mode) than for rotations about the tetragonal axis (Ai g mode), in agreement with experimental 
results. The FE mode, on the other hand, is found to be softer when polarized perpendicular to 
the tetragonal axis (E u mode) than parallel to it (A2 U mode). The sign of this frequency splitting 
is consistent with the experimentally reported anisotropy of the dielectric susceptibility and other 
evidence. Finally, we present a discussion of the influence of various types of structural distortions 
on the FE instability and its anisotropy. 

PACS numbers: 77.84.Dy, 77.80.Bh, 77.80.-e, 63.20.-e 
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I. INTRODUCTION 

First-principles calculations are proving to be one of 
the most powerful tools for carrying out theoretical stud- 
ies of the electronic and structural properties of ma- 
terials. A particularly successful application of this 
technique has been its use in understanding the per- 
ovskite ferroelectric compounds. These materials have 
important technological applications because of their 
switchable macroscopic polarization and their piezoelec- 
tric properties. They are also attractive objects of fun- 
damental study because of the rich variety of phase dia- 
grams that they display as a function of temperature. At 
high temperature, the ABO3 perovskites retain full cubic 
symmetry. However various structural phase transitions 
take place as the temperature is reduced.u For example, 
BaTi03 and KNbC>3 undergo phase transitions from the 
cubic paraelectric (PE) phase to a succession of tetrago- 
nal, orthorhombic, and finally rhombohedral ferroelectric 
(FE) phases. In contrast, PbTiC-3 displays only a single 
transition, from the cubic PE phase to a tetragonal FE 
phase. In NaNbC>3 and PbZrC>3, non-polar antiferrodis- 
tortive (AFD) or antiferroelectric (AFE) transitions take 
place, associated with different types of tilts of the oxy- 
gen octahedra, in addition to the FE transitions. It is 
understood that the two types of transitions result from 
the condensation of soft phonon modes at the Brillouin 
zone boundary with q and at the zone center with 
q = 0(Ref.|). 

The low-temperature behavior of SrTiC>3 has been an 
attractive subject for experimental and theoretical study. 
SrTiC>3 behaves as an incipient ferroelectricu (similar to 
KTaOs) in the sense that it has a very large static dielec- 



tric response and is only barely stabilized against the 1 
densation of the FE soft mode at low temperature.! 
As the temperature is reduced, SrTiC>3 first undergoes 
a transition from the cubic to a tetragonal AFD phase 
at 105Kj3 but this transition is of non-polar character 
and has little influence on the dielectric properties. The 
static dielectric response closely obeys a Curie- Weiss law 
of the form of e ~ (T — T c ) _1 at temperatures above 
about 50K, but the divergence at a critical temperature 
T c ~ 36K thaL would be expected from this formula is 
not observedaM Instead, the susceptibility saturates at 
an enormous value of ~ 2 x 10 4 as T approaches zero. 
Because the system is so close to a ferroelectric state, it 
is not surprising to find that it can be induced to be- 
come ferroelectric, either by the application of electric 
fielcB, uniaxial stressH or by the substitution of Ca ions 
on the Sr sublattice.BI Finally, thq-StliOs system also 
displays puzzling phonon anomalieai3~li3 and electrostric- 
tive responseEj in the low-temperature regime. 

This peculiar behavior, especially the failure of the 
system to condense into a FE phase at low T, has 
been the subject of considerable theoretical study and 
speculation.EjLj Recently, efforts have focused on the so- 
called "quantarn paraelectric state" postulated by Muller 
and BurkardJiH who suggested that quantum fluctuations 
of the atomic positions could suppress thepEE-transition 
and lead to a stabilized paraelectric state.tlrEJ This hy- 
pothesis has received dramatic support from a recent ex- 
periment showing that isotopically exchanged SrTi 18 03 
appears to become ferroelectric at 23K,E3 suggesting that 
normal SrTi 16 03 must be very close indeed to the ferro- 
electric threshold. 

First-principles calculations have already contributed 



1 



significantly to the ujftderstaxiding of the structural 
properjtifls-.^ SrTi0 3 .BBm23 Calculations of two 
groupsEilcJ'CJ confirmed that SrTiOa, in its high- 
symmetry cubic structure at T = 0, is unstable to both 
FE and AFD distortions when the atomic coordinates 
are treated classically. Using classical Monte Carlo sim- 
ulations on an effective-HamiltoniarES fitted la the first- 
principles calculations, Zhong and Vanderbiltcil predicted 
that SrTiOa would first undergo the AFD transition at 
about 130 K, and then a further transition into a state 
with simultaneous AFD and FE character at 70 K. An- 
harmonic interactions between the AFD and FE modes 
were found to be competitive, in the sense that the pres- 
ence of the AFD distortion was found to reduce the FE 
transition temperature by about 20%. The same authors 
later showed that, when a quantum-mechanical treat- 
ment of the atomic positions was included via a quantum 
path-integral Monte Carlo simulation, the AFD transi- 
tion temperature was shifted very close to the experimen- 
tal one at 105 K, and the FE transition was suppressed 
down tcUjhe lowest temperatures that could be studied 
(~5 K),El consistent with the experimental absence of 
a transition. On the other hand, LaSota and coworkers 
have recently performed a first-principles calculation of 
the ground state structural properties and interactions 
between the FE and. AFD instabilities in SrTiOs using 
an LAPW approach .E3EZI These authors found that the 
AFD tetragonal structure is stable against the FE dis- 
tortions, indicating no interaction between the AFD and 
FE modes. |— The conclusion is in contrast to the previ- 
ous theory.Ell However, LaSota et al. made certain ap- 
proximations to the FE eigenmodes, as will be discussed 
below. 

An additional motivation for a detailed theoretical 
study of SrTiOa is the opportunity to make contact with 
the remarkably systematic experimental study of Uwe 
and Sakudo.tj These authors made careful measurements 
of the anisotropic dielectric susceptibilities and Raman 
mode frequencies as a function of uniaxial stress applied 
along different crystal orientations. They also fitted their 
results, plus those of previous experimental studies, to 
obtain a phenomenological description of the couplings 
between the AFD, FE, and strain degrees of freedom of 
the crystal. In particular, their fit contains an anisotropic 
coupling which, in the tetragonal AFD phase, tends to 
favor FE distortions that are perpendicular to the tetrag- 
onal axis over those that are parallel. However, previ- 
ous theoretical work has given an unclear picture of this, 
anisotropy. On the one hand, Vanderbilt and ZhongEa 
found that the interaction between the FE and AFD 
modes, which was mainly through the on-site anharmonic 
coupling, would tend to favor FE modes polarized per- 
pendicular to the AFD tetragonal axis, in accord with 
experiment. On the other hand, the Monte Carlo calcu- 
lations previously referred toEil indicated a sequence of 
transitions with decreasing temperature in which the FE 
order parameter of a z-polarized mode was found to de- 
velop before the x or y-polarized ones, indicating that 



the z-polarized mode-jsoes soft first. Finally, the recent 
work of LaSota et aZ.oEO suggests that there is very little 
anisotropic coupling at all. In view of these apparently 
conflicting theoretical results, we felt it worthwhile to 
clarify this situation by carefully studying the anisotropy 
of the FE distortion energy in the AFD ground state. 

With these motivations, we have carried out a thor- 
ough analysis of the ground-state structural and dynam- 
ical properties of tetragonal SrTiOs using first-principles 
density-functional calculations. This approach is based 
on a classical treatment of the nuclear motions, and so 
is obviously unable to take into account the quantum 
zero-point motion of the ionic positions which becomes 
critically important at low temperature. Nevertheless, 
from such a calculation one is still able to compute in- 
teraction parameters for comparison with experiment, to 
identify the effects which tend to suppress the FE insta- 
bility in the presence of the AFD state, and to obtain 
qualitatively a picture of the dielectric anisotropy con- 
nected with the splitting of the differently polarized FE 
modes in the AFD state. 

The calculations are carried out using a plane-wave 
basis and ultrasoft pseudopotentialsBa The theoretical 
equilibrium AFD structure is obtained by minimizing the 
energy with respect to cell volume, c/a ratio, and in- 
ternal parameters. Frozcn-phonon calculations are then 
used to obtain the frequencies of T-point and i?-point 
phonon modes, including FE soft modes, in the AFD 
ground state. For this purpose, we make use of a point- 
group symmetry analysis to reduce the complexity of 
the distortions that need to be studied. In order to in- 
terpret the results in terms of a phenomenological de- 
scription involving FE and AFD mode distortions and 
strains, we use an approach similar to that underlying 
the effective-Haffliltonian scheme first developed-for the 
BaTi0 3 systemEa and later applied to SrTi0 3 x9 That 
is, we use the LDA calculations to compute the values 
of the Taylor expansion coefficients of the total energy 
with respect to these distortions, and compare with the 
experimental determinations of Uwe and SakudoB 

The rest of the manuscript is organized as follows. In 
Sec. II we briefly describe the technique employed for 
the first-principles calculations. In Sec. Ill we present 
and discuss the results of the calculations. We begin 
with the determination of the theoretical tetragonal AFD 
structure, and then proceed to study the energies of AFD 
and FE distortions about this reference structure, with 
special attention to the anisotropics of the AFD and FE 
mode frequencies. Finally, we conclude in Sec. IV. 

II. THEORETICAL DETAILS 

Our ab-initio plane-wave pseudopotential calcula- 
tions are based on the Hohenberg-Kohn-Sham density- 
functional theory (DFT) within the local-density approjx. 
imation (LDA). Ultrasoft Vanderbilt pseudopotentialsEa 
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FIG. 1. Sketch of the TiC>2 layer in the unit cell of the 
tetragonal AFD structure, illustrating rotation of octahedra 
about the z axis, x-y and x'-y' are the coordinate frames in 
the tetragonal cell and the original cubic cell, respectively. 

are used, with the 0(2s), 0(2p), Ti(3s), Ti(3p), 
Ti(3d), Ti(4s), Sr(4s), Sr(4p), and Sr(5s) states in- 
cluded in the valence. The exchange-correlation en- 
ergy is of the Cenedey- Alder form with Perdew-Zunger 
parameterization.LJLJ A conjugate-gradient minimiza- 
tion schemea is used to minimize the Kohn-Sham en- 
ergy, using a plane-wave cutoff of 30 Ry for all calcula- 
tions. Unless otherwise stated, our calculations are car- 
ried out at [the theoretical equilibrium lattice constant 
of 7.303 a. u.Ej which is ~1% less than the experimental 
value of 7.365 a. u., the discrepancy representing the in- 
herent LDA error. 

Cubic SrTiOa has a simple cubic 5-atom unit cell with 
a common lattice parameter a along the [100], [010], and 
[001] directions. We will briefly discuss some calculations 
carried out for a doubled unit cell corresponding to the 
condensation of a soft AFD mode at the (110)7r/a or 
M point of the Brillouin zone (BZ) boundary. However, 
most of our attention will be focused on the ground-state 
tetragonal phase obtained by freezing in an AFD phonon 
mode at the (lll)7r/a or R point of the BZ boundary. 
This triply-degenerate phonon mode corresponds to the 
rotation of the TiO@ octahedra in opposite directions 
from one cubic unit cell to the next, followed by a small 
tetragonal strain. (Note that it is conventional to label 
the phonon modes with respect to the simple-cubic BZ, 
even when they condense to lower the symmetry.) Taking 
the rotation to be about the z axis, we adopt a 10-atom 
tetragonal unit cell with lattice vectors of length >/2a, 
i/2a, and c along the [110], [110], and [001] directions, 
respectively. (That is, in our convention, c/a is close to 
1, not 1/V2-) The rotation of the oxygen atoms in the 
Ti— O plane is shown in Fig. [j]. Throughout this paper, 
we will use x' and y' to denote the original cubic direc- 
tions ([100] and [010], respectively), while x and y are 
taken as parallel to the tetragonal lattice vectors along 
[110] and [110], respectively. That is, the x-y frame is 
rotated by 45° relative to the x'-y' frame (z axes are 
congruent). 

In all cases, we use a k-point set that is equivalent to 
the 6x6x6 Monkhorst-pack mesl£3 in the BZ of the sim- 
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FIG. 2. Calculated values (symbols) and fits (curves) of 
the total energy per 10-atom cell as a function of octahedral 
rotation angle 9 Z , computed at the frozen theoretical cubic 
lattice constant a=7.303a.u., for both J?-point (diamonds) 
and M -point (circles) octahedral rotation modes. 

pie cubic cell, corresponding to 108 k-points in the full 
BZ of the tetragonal cell. The irreducible BZ then con- 
tains 6 k-points for the undistorted cubic structure; 10 
k-points for the tetragonal ground-state structure, with 
or without additional A\ g or A^ u mode displacements; 
and 20 k-points for the tetragonal structure with addi- 
tional E u mode displacement. 

III. RESULTS AND DISCUSSIONS 

A. AFD instability in cubic unit cell 

To establish notation, we let the energy of an AFD 
phonon mode per 10-atom cell of the cubic perovskite 
structure be expanded up to fourth order in <j) z , 

E = E + ^ 2 Z + A x( j ) j , (1) 

where ip z = (a/2) sin 0^ is the magnitude of the oxygen- 
atom displacement associated with the rotation of the 
oxygen octahedra around the [001] axis, as shown in 
Fig. [l]. In Fig. || we show the computed values of the 
total energy versus rotation angle for AFD modes at 
both the M and i?-points of the BZ (corresponding to 
in-phase or out-of-phase rotations in neighboring planes 
of octahedra along z respectively). These were com- 
puted at zero strain, i.e., with the lattice vectors fixed 
to be those of the theoretical equilibrium cubic structure 
(a=7.303a.u.). As can be seen, the computed total en- 
ergy versus rotation angle can be fitted very well by the 
quartic Eq. ([j]). Defining the mode stiffness n=d 2 E /dcj) 2 zl 
we find that k < for both M- and i?-point modes as 
shown in Fig. |^. Nevertheless, the magnitude of k for the 
M -point mode is only ~10% that of the i?-point mode, 
indicating that the instability at the R point is much 
stronger than that at the M point. Consequently, for 
the remainder of this paper we will limit our discussion 
to i?-point distortions only. 
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As can be seen from Figs. ^]and 3, the equilibrium octa- 
hedral rotation angle is found to be (9 Z =5.5°, significantly 
larget-than the zero-temperature experimental value of 
2.1°I3 Since the theoretical equilibrium lattice constant 
(7.303 a. u.) is somewhat smaller than the experimen- 
tal one, we also carried out similar total-energy calcu- 
lations at the extrapolated zero-temperature (7.365 a. u.) 
and roomptemperature (7.38 a.u.) experimental lattice 
constants.Ea The results shown in Figure 3 confirm that 
increasing the lattice constant or crystal volume tends 
to suppresS|-t|he, AFD instability, as expected from pre- 
vious work.ESE2l However, the resulting variation of the 
equilibrium rotational angle is too small to explain the 
experimental observation, changing only marginally to 
4.89° and 4.69° at the zero- and room-temperature ex- 
perimental lattice constants, respectively. 

These results demonstrate that the underestimate by 
~ 1% of the lattice constant by the LDA is not the pri- 
mary factor responsible for the theoretical overestimate 
of the rotation angle. Moreover, we shall see in the next 
subsection that the inclusion of strain relaxation effects 
only acts to increase (slightly) the theoretical equilibrium 
rotation angle. Thus, we think that the smaller observed 
value of the AFD rotation angle can most likely be at- 
tributed to the quantum fluctuations associated with the 
motion of the oxygen atoms. This effect is not included 
in the theory, and should act to reduce the amplitude 
of symmetry-breaking distortions. While previous work 
has indicated that the quantum fluctuations should have 
a weaker effect on the AFD modes than upon the FE 
ones,E2l the effect on the AFD modes could still be quite 
significant. An alternate possibility is simply that the 
underestimate is a result of LDA error not associated 
with the lattice constant. In any case, we have chosen 
to complete our theoretical investigations by consider- 
ing distortions about our theoretical ground-state AFD 
structure, keeping in mind that the results should be in- 
terpreted with the overestimate of the rotation angle in 
mind. 



B. AFD modes in the tetragonal structure 

To study the ground-state tetragonal structure, the 
lattice strains also need to be taken into account. We 
adopt the usual Voigt notation Xi for the strain tensor, 
but set X4 = X5 = xq = because such off-diagonal shear 
strains will not enter into our considerations. We take 
i=l, 2, 3 corresponding to the x', y', and z pseudocubic 
axes for both strains Xi and rotations fa . Expanding the 
energy up to quartic order in AFD amplitudes, quadratic 
order in strain, and leading order in the strain- AFD cou- 
pling, the symmetry-allowed contributions are 
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FIG. 3. Calculated values (symbols) and fits (curves) 
of the force experienced by an oxygen atom as a func- 
tion of octahedral rotation angle (9 Z , computed at the 
theoretical cubic (a=7.303 a.u., diamonds), the experi- 
mental T=0 (a=7.365 a.u., circles), and the experimental 
room-temperature (a=7.38a.u., squares) cubic lattice con- 
stants. Arrows indicate the theoretical and experimental equi- 
librium rotation angles of 5.5° and 2°, respectively. 



-c n x i + C 12 XiXj 
i i<j 

hi Y X i ( t ) 'i ~ ^12 Y Xi< ^T' 



(2) 



i<3 



We choose the tetragonal ground state to be oriented 
along the z-axis, with x\ = X2 7^ 0, fa — fa — 0, and 
fa =/= 0. In this ground state, phonon modes correspond- 
ing to additional oscillations of fa belong either to the E g 
(i=l or 2) or the Ai g (i=3) representation of the tetrag- 
onal Z?4/j point group. For symmetry-preserving (Ai g ) 
distortions, it is convenient to re-express the three strain 
components in terms of a volume strain x and a shear 
strain v according to 



X\ = X2 — x — v , 



X3 = x + 2v 



(3) 



In terms of these variables, the second and third lines of 
Eq. (0) can then be rewritten as 



£•^(2,0) = -oaT + 



(4) 



and 



^coupling^ ^ .) = ^-^2 fal - lv (2fai -fa{- fal) , 

i 

(5) 



where 



a = c\\ + 2ci2 , 
/3 = en - C12 , 

T) = fon + 26i2 , 

7 = 6n - 612 . 



(0) 
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TABLE I. Parameters in the effective Hamiltonian of 
SrTiOa. See text. 
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To find the equilibrium rotation angle in the A\ g 
ground state, we hold 4>i—4>2—^ and minimize E in 
Eq. (^) with respect to x and v at a fixed (f>3 — <p z - The 
minimizing values are 

3/r z 



(7) 



,,cq 



Substituting into Eq. (g 

s) = Eq 



E 



where 



.4 



x 



— nq 

2 V 



rf_ 
6a 



A 



11 
3/3 



(8) 



(9) 



Thus, when the strain relaxation is taken into account, 
the equilibrium rotation angle is given by 



4A 



x 



(10) 



We determine all the interaction parameters k, A Xl 
A™, en, C12, &n, and 612 via a series of finite-difference 
calculations of total energies and forces within the LDA. 
Table I lists our results and compares them with the cor- 
responding values determined by Uwe and Sakudo by 
fitting to experiments (all units have been converted jta 
atomic units). We find very good agreement overall.E3 
The fact that Ax < A x implies that the inclusion of 
strain relaxations strengthens the AFD instabilities at 
anharmonic order. The equilibrium rotation angle in- 
creases to 6.0° (compared with 5.5° for the cubic strain 
state), while the equilibrium values of x and v are found 
to be -0.10% and 0.23%, respectively. 

In the tetragonal AFD ground state, the frequencies of 
the soft phonon modes associated with additional rota- 
tions of the oxygen octahedra are given by evaluating 



d 2 E 



(11) 



where — 4mo is the mass factor associated with the 
oxygen rotational mode and the derivative is to be eval- 
uated under conditions of fixed strain at the equilibrium 
structure (i.e., at the equilibrium values of x, v, and (f> z ). 
Eq. |ll| gives the frequencies of the E g and Ai g modes to 
be, respectively, 



uj\ = -K(A x /2m^A x ) , 
-2KA x /m^A x , 



(12) 



where A x 



A™ + Y/[3. The values of the two fre- 
quencies are 45 cm -1 and 130.7cm -1 respectively, so that 
the frequencies of the A\ g and E g modes are roughly in 
the ratio 3:1. This is consistent with observed ratios of 
^2.5:1 in pressure-dependent experipapntsEJ and ~3:1 in 
temperature-dependent experimentsE3. 

Up to this point, the analysis has been done at the the- 
oretical equilibrium lattice constant. However, it is well 
known that the LDA tends to underestimate the lattice 
constants of perovskites by ~l%.c3 Moreover, past expe- 
rience has shown that the displacement patterns associ- 
ated with soft modes ma y de pend critically on the lat- 
tice constant and strains .E§E3 To take these effects into 
account, we adopted a strategy of applying a negative 
hydrostatic pressure to .the lattice to restore the experi- 
mental lattice constant .lj Using Eq. (|^) and minimizing 
the Gibbs free energy 



G = E + 3xP 
with respect to x at fixed pressure P, we find 



p 

a 



(13) 



(14) 



while v is independent of pressure. Then 



G(<^,P) = ~K eff 



3_P 2 
~2a 



-Z- (15) 



where the effective harmonic coefficient is 

„ P 
K c ft = k + 2rj — . 

a 

Thus, at fixed pressure, the equilibrium rotation angle is 



Kcfi 



4A 



X 



(16) 



As one can see, the harmonic coefficient K e g depends 
upon the external pressure variable. It would thus be 
possible, in principle, to adjust P so as to fit the re- 
sulting rotation angle to the experimental angle of 2.1°. 
However, the pressure needed to achieve this, — 14.4GPa, 
would expand the lattice constant to 7.48 a. u., which 
is much larger than the experimental value. Instead, 
we adjust P so as to fit the experimental lattice con- 
stant. That is, we adjust P so that the volume strain is 



x = (a CX p-atheo)/atheo = 0.849%, where a, 



cxp~ 



=7.365 a. u. 
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and a t heo are the zero-temperature experimental and the- 
oretical lattice constants, respectively. Substituting (|l6| ) 
into (|l4|), we obtain 



TABLE II. Character table for point group D4 



x(P) = x 



P 

OieS 



(17) 



where 



x = 



KTj 



YlaAx ' 



a e ff = a 



1 



6aA 



x 



(18) 



Inserting x = 0.849% leads to P = -5.26GPa. The 
strains along the tetragonal and planar axes are found to 
be 1.135% and 0.705%, respectively. Relative to a cubic 
cell with the experimental lattice parameter 7.365 a. u., 
the tetragonal cell we obtained is thus expanded along 
[001] while compressed along the [100] and [010] direc- 
tions. In this circumstance, the equilibrium rotation an- 
gle is found to be 4.93°. 

Under these conditions, the A\ g and E g soft-mode fre- 
quencies of Eq. ([l2]) now become 



ujf = u\ = -K cS (A x /2m 4> A x ) 
ujj = -2K e fiA x /m<j,Ax ■ 



(19) 



From the coefficients in Table A x /Ax = lA87 and 
2A X /A X =1.11, so that clearly \uj(E g )\ < \uj(A lg )\. 
In fact, the two frequencies are calculated to be 
uj (E g )— 37 cm^ 1 and u>(Axg)=109 cm' 1 . (For compari- 
son, w(Aig)=124cm~ 1 in Ref. ^7], while the measured 
oj(E g ) and u(A%g) in Ref. || are 15 cm -1 and 48 cm -1 
respectively) . 

To summarize the results so far, we have found that 
the AFD mode condenses at the R point of the cubic BZ, 
associated with a triply-degenerate phonon of sym- 
metry. As a consequence of the transition from the cubic 
to the tetragonal state, the degenerate i?-point modes 
split strongly into an A\ g singlet and an E g doublet, the 
latter having a softer frequency than the formei^. This is 
in good qualitative agreement with experiment u 

The expansion approach used above for the AFD 
modes makes an implicit assumption that the phonon 
eigenvectors from the cubic structure are a good approx- 
imation to those in the tetragonal AFD structure. For 
the AFD modes, where the anisotropy is large, we do not 
expect this approximation to be at all serious. However, 
our next task will be to analyze the FE mode anisotropy 
in the AFD state. As will be seen below, this turns 
out to be much more delicate than for the AFD modes. 
Thus, we have chosen to take a different approach for 
the FE modes, in which the normal modes in the AFD 
ground state are directly computed. The symmetry anal- 
ysis needed to do this is given in the next subsection, and 
the FE mode analysis is then given in the concluding sub- 
sections. 
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— 1 


B 2 g 


1 


—1 


1 


—1 


1 


—1 


—1 


1 


—1 


1 


Eg 


2 





-2 





2 








-2 








A lu 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


A 2u 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


1 


1 


Bin 


1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


B 2u 


1 


-1 


1 


1 


-1 


1 


1 


-1 


1 


-1 


E u 


2 





-2 








-2 





2 









C. Symmetry analysis of normal modes in AFD 
tetragonal phase 

In this section, we present some details of the point- 
group symmetry analysis of the normal modes in the 
AFD tetragonal structure, needed for the calculation of 
the frequencies of transverse optical (FE) modes at the 
Brillouin zone center. 

To harmonic order, the displacement energy can be 
expressed as 



E = l 
2 



(20) 



where uf is the displacement of sublattice i in Cartesian 
direction a, and the force constant matrix $ obeys the 



symmetry conditions $ 



a, {3 
i,3 



cDff and V -!> , 



0. The 



dynamical matrix is just related to the force constant 
matrix by a diagonal mass tensor. 

It is well known that the vibrational modes at a 
given k-point in the BZ of a crystal transform accord- 
ing to the corresponding irreducible representations of 
the symmetry group for that k-point. Such an analysis, 
which has previously been used to construct the force- 
constant matrices-|for the FE modes in the cubic per- 
ovskite structureE3 and foi all modes in the tetragonal 
FE structure of PbTiOsJfj is applied here to the zone- 
center modes in the AFD tetragonal structure of SrTiC>3. 
When SrTiC>3 condenses from the cubic into the tetrago- 
nal AFD phase, the point group lowers from Oh to D^h- 
Because we are interested in zone-center modes, the sym- 
metry group of k is just the D^h point group itself, which 
contains 16 symmetry operations that can all be gener- 
ated from a fourfold rotation C4 and two mirror reflec- 
tions ah and ad- The character table is shown in Table 
II. There are 10 irreducible representations (irreps), of 
which two are two-dimensional. The AFD soft modes 
originating from the cubic T2s(i?) phonons now belong 
either to the A\ g or E g irreps, depending on whether the 
octahedron rotation axis is parallel or perpendicular to 
the tetragonal axis, respectively. Similarly, the modes 
originating from cubic Ti5 FE modes are now either A2 U 
or E u , depending on whether the polarization is parallel 
or perpendicular to the AFD axis. 
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TABLE III. Symmetry analysis of the normal modes in 
tetragonal AFD structure. 



irreps 



dim. 



Sr 



Ti 



Oi- 



o. 



A 2g 

Big 
B 2g 

E 9 

Aiu 

A 2u 

Biu 

B 2u 

E u 



SzR 
S(x,y)R 
SzF 

S(x,y)T 



TzR 
TzY 



T(x,y)(R) 



i? a 
R a 
R a 
i? a 
{E,0)zR 

EzT 

OzT 
E{x,y)T 
0(x,y)T 



SzR 

3(x,y)R 
3zF 

3(x,y)T 



'From decomposition of (E,0)(x,y)R. 



To label all of the displacement patterns associated 
with the 10-atom cell (30 degrees of freedom), we use 
the notation TaK to denote a displacement associated 
with atom type T in pseudocubic Cartesian direction a 
and having "phase relation" K. The five atoms types are 
abbreviated as t S' for Sr, £ T" for Ti, '3' for oxygen atoms 
making Ti-0 chains in the z direction, and '1' and '2' for 
oxygen atoms in TiC>2 x-y planes. The "phase relation" 
K is either T' or 'i?' depending on whether the two atoms 
of the same type in the 10-atom cell move in-phase or out- 
of-phase (that is, whether they originate from T-point 
or i?-point modes of the parent cubic structure). Note, 
however, that some individual displacements contribute 
to more than one irrep (e.g., lzT and 2zT contribute to 
both the Am and Biu irreps). Thus, for the two in-plane 
oxygens we introduce alternative "type" designations E 
('even') and O ('odd') in place of '1' and '2', where 



1 / \ 

U-EaK — —^{UlaK + U2aK) 



UQaK — —^(u2aK ~ U\ a K) 



(21) 



(six degrees of freedom for each atom type), classify- 
ing them according to the irreducible representations to 
which they belong. Note that the four (E,0)(x,y)R 
modes do not have a simple one-to-one correspondence 
with the Ai g , Ai g , R\ g , and Bi g modes to which they 
give rise; these are indicated in the table with just the 
notation R. 

Since, we have access to the Hellmann-Feynman 
forcescil in our first-principles ultrasoft-pseudopotential 
approach, it is convenient to compute the force-constant 
matrix elements from finite differences as 



a,0 



Ff_ 



(22) 



where F is the force that results from a sufficiently small 
displacement u. We can use the symmetry analysis to 



TABLE IV. Zone-center 


transverse 


optical 


phonon fre- 


quencies in the cubic structure. (All units in cm 




rpoint 


TOl 


T02 


Tf>Q 

1 Uo 


Current (7.303 a.u.) 


42 


168 


o4y 


Current (7.365 a.u.) 


94i 


151 


521 


PW a (7.30 a.u.) 


41i 


165 


546 


LR b (7.412 a.u.) 


lOOi 


151 


522 


Expt. c 90K 


42 


175 


545 



a Ref. W2: Plane- wave pseudopotential method. 

b Ref. yA: LAPW linear response method. 

c Ref. [S Fitted from experimental infrared reflection spectra. 



identify the set of sublattice displacements that may par- 
ticipate in a given normal mode, and to calculate the 
forces that arise at first order with each such displace- 
ment. For example, for the A^ u FE mode, we find the 
four displacements of types SzT, TzT, EzT, and 3zr may 
participate. For each, a displacement amplitude of 0.2% 
of the lattice constant is chosen so that the harmonic 
approximation is still well satisfied. From each such cal- 
culation, the resulting force vector is projected onto the 
same set of four displacements, thus building up the 4x4 
force-constant matrix. This matrix is then symmetrized 
and diagonalized. In a similar way, the FE E u mode 
is represented in a 6x6 subspace with basis SxT, TxT, 
TyR, SxT, ExT, and OyT. 



D. FE instability in the AFD tetragonal phase 

Previous work of Zhong and VanderbiltS and LaSota 
et alE-l has indicated that cubic SrTi03, in the absence 
of any AFD distortion, shows a FE instability. We have 
confirmed this result here by performing frozen-phonon 
calculations and building up the 4x4 force-constant ma- 
trix for Ti5 modes polarized along z (essentially the same 
procedure outlined for the A2 U modes at the end of the 
last subsection, except performed for the cubic 5-atom 
cell). The non-zero eigenvalues of the corresponding dy- 
namical matrix are given in Table IV for both the theoret- 
ical equilibrium (7.303 a.u.) and expanded experimental 
(7.365 a.u.) lattice constants, together with other the- 
oretical and experimental results for comparison. Our 
results at the theoretical volume indicate that the FE 
soft phonon mode frequency remains real, although it 
was imaginary according to the earlier linear-response 
and plane-wave calculations which used slightly different 
lattice constants.E-rc3 However, at the zero-temperature 
experimental volume, the frequency is found to be imag- 
inary, indicating that the FE instability is indeed very 
sensitive to the crystal volume. Increasing the crystal 
volume enhances the FE instabilities, in agreement with 
experiment Qcil This is in contrast to the AFD instabil- 
ity, which weakens with increasing volume as shown in 
Sec. III. A. These results suggests that there is an inher- 
ent competition between AFD and FE modes. 
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TABLE V. Calculated eigenfrequencies and displacement 
patterns for the A2« mode in tetragonal SrTi03. Columns 
are labeled by mode eigenfrequency; imaginary frequency in- 
dicates an unstable mode. 



90i cm 



157 cm 



515 cm 



SzT 
TzV 
EzY 
?>zT 



-0.172 
-0.569 
0.636 
0.492 



-0.700 
0.641 
0.285 
0.128 



-0.036 
0.058 

-0.583 
0.809 



TABLE VI. Calculated eigenfrequencies and displacement 
patterns for the E„ mode in tetragonal SrTiOa. Columns are 
labeled by mode eigenfrequency. 



96i cm 



240 cm" 



419 cm" 



SxT 
TxF 
TyR 
ExF 
OyV 
3aT 



0.217 
0.524 
0.086 
-0.668 
-0.034 
-0.472 



0.028 
-0.015 

0.014 
-0.518 

0.497 

0.694 



0.013 
-0.077 

0.995 

0.055 
-0.008 

0.024 



To understand the behavior of the FE modes in the 
zero-temperature AFD tetragonal structure, we next per- 
form frozen-phonon calculations for this structure. We 
use the AFD state described at the end of Sec. III.B, 
i.e., P = -5.26GPa, V =798.989a.u. 3 , c/a = 1.004, and 
</> 2 =4.93°. Both E u and Am modes are calculated us- 
ing the method presented in Sec. III.C. Tables V and 
VI present the resulting eigenfrequencies and eigenvec- 
tors for the three lowest modes of each symmetry. It can 
be seen that both the E u and A2 U soft-mode frequencies 
are imaginary (96i and 90icm _1 , respectively). Total- 
energy calculations are then performed with the corre- 
sponding eigenmodes at different mode amplitudes. The 
double-well energy curves are as shown in Fig. ^. One can 
see that the well depth of the E u mode is substantially 
greater than that of the A2 U mode, indicative of a much 
stronger instability for a distortion of E u symmetry. 

With respect to the FE mode frequencies, the com- 
parison with experiment is problematic because the the- 
oretical values are imaginary (indicating an instability) 
while the experimental values are not (consistent with a 
non-FE T=Q ground state) . As indicated in the Introduc- 
tion, the experimental stabilization of the non-FE AFD 
structure is understood to result from quantum fluctua- 
tions of the atomic coordinates. It should be emphasized 
that the purpose of the present calculations is to study 
the low-temperature structural instabilities in a classi- 
cal framework, i.e., in the absence of quantum fluctua- 
tions. Thus, the existence of the FE instabilities in our 
DFT-LDA ground state calculations, takeii-together with 
earlier quantum Monte Carlo simulations tend to cor- 
roborate the hypothesis that the quantum fluctuations 
are responsible for the experimental absence of a real FE 
phase. 
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FIG. 4. Total energy vs. the FE distortion amplitude at 
AFD tetragonal phase. Circles and solid curve are for E u 
mode; diamonds and dashed curve are for A% u mode. 



Nevertheless, we can make the following comparisons. 
In the stress-induced Raman scattering measurement of 
Uwe and Sakudo, the E u and Am soft mode frequencies 
were found to be 9.1±0.6cm _1 and 19±lcm _1 respec- 
tively (this appreciable splitting between the two modes 
has frequently been overlooked) .□ The two squared fre- 
quencies produce a difference of 3x 10 2 cm" 2 . Our calcu- 
lation shows the difference between the two frequencies 
to be 6 cm -1 , roughly the same order of magnitude as the 
difference found experimentally. However, the theoreti- 
cal difference of squared frequencies is 10xl0 2 cm~ 2 , or 
roughly three times larger than the experimental value. 
This is simply because of the large magnitude of the 
imaginary frequencies. 

In summary, we found imaginary frequencies for both 
the FE E u and A2 U soft modes in the AFD tetragonal 
phase. There is an apparent splitting between the two 
modes, ui 2 (E u ) < ui 2 (A2 U ), suggesting that the FE struc- 
ture of E u symmetry is more energetically favorable than 
the A2 U one. This result is consistent with the fact that 



the A2 U FE mode i, 
tcring experiments! 
that of the E u phonon.B 



ess easily observed in neutron scat- 
bgcause its energy is higher than 



E. Influences of structural distortions on the 
stabilities 

In the previous subsection, we have calculated the FE 
phonon frequencies for either a cubic structure in the ab- 
sence of the AFD distortion, or for a tetragonal AFD 
phase as observed experimentally. To understand the 
interaction of different distortions and their roles in af- 
fecting the FE instabilities, we performed frozen-phonon 
calculations for the cubic reference structure at the ex- 
perimental lattice constant both with and without the 
tetragonal strain, and with and without the AFD rota- 
tion. In Tabic VII we present the results for the FE 
phonon frequencies for each of these scenarios. 

Observe that the two FE modes instabilities with dif- 
ferent symmetries depend on the strain distortion and 
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TABLE VII. Calculated E u and phonon frequencies at 
different distorted structures for SrTiOa. x, v and 9 Z denote 
uniform volume strain, shear strain, and rotation angle about 
the z axis, respectively. 



Distortions 



FE symmetry 



X 


V 


e z 


A 2 u 


E u 











42 


42 


0.00849 








94* 


93* 


0.00849 


0.00143 





112* 


83* 


0.00849 


0.00143 


4.9° 


90*' 


96* 


0.00849 





4.9° 


63*' 


107* 



AFD distortion in an opposite sense. For the A 2u mode, 
a non-zero shear strain increases the FE instability while 
an AFD rotation distortion reduces it. The E u mode, 
however, depends on the two distortions in the opposite 
way. Thus, the final sign of the frequency splitting (i.e., 
the sign of the anisotropy) ultimately depends sensitively 
on a delicate partial cancellation of the two contribu- 
tions. The final result is that the E u mode is energeti- 
cally slightly more unstable than the A 2u mode. 

Finally, we have tested whether the use of the cubic 
eigen-modc distortion in place of the true tetragonal one 
is a good approximation for understanding the coupling 
of the FE and AFD instabilities. We calculated the ex- 
pectation value (£, c \&t\£,c) , where $t is the dynamical 
matrix in th e tet ragonal ground state structure (fourth 
row of Table VII ) and £ c is the cubic eigen-mode distor- 
tion vector. The resulting imaginary frequencies are 85.4* 
and 86.9* cm -1 for the A 2u and E u modes respecti vely , 
compared with values of 90* and 96icm~ 1 in Table VII. 



Thus, the splitting is much smaller when the expectation 
value is used. This makes it clear that the nature of the 
FE anisotropy is quite sensitive to the tetragonal lattice 
strain associated with the AFD rotation distortion. This 
suggests a probable explanation for the fact that the clas- 
sical MC simulations of Ref . 20 , which are based on use 
of a cubic mode eigenvector, predicted the anisotropy in- 
correctly - i.e., the A 2u mode was found to go soft before 
the E u one. 



IV. CONCLUSION 

In this work, we have investigated both the AFD 
and FE structural instabilities in the tetragonal phase 
of SrTiC"3. A unique aspect of this work is that we 
have studied the FE frequencies carefully using the ex- 
act eigenmode distortion obtained from the ground-state 
tetragonal structure, whereas previous studies have made 
the approximation of using the cubic eigenmode distor- 
tion instead. We show that the instabilities have a sen- 
sitive dependence on the crystal volume. We also found 
that the existence of the FE instabilities are affected by 
the coupling to the shear strain and the rotational dis- 
tortion. Both types of distortion contribute, but with 



opposite sign, so that it is a subtle cancellation between 
them that determines the splitting of frequencies between 
the E M and A 2u modes. The E u mode is found to be the 
more unstable of the two in the ground-state AFD phase. 

The degree of tilt of the oxygen octahedra is still over- 
estimated in our calculations, compared to the experi- 
mental results. We attribute this to the quantum fluc- 
tuations of the oxygen coordinates, which tend to sup- 
press the average AFD distortion amplitude in the ex- 
periments. The AFD and FE instabilities show opposite 
trends with an increase of crystal volume (weakening and 
strengthening, respectively). However, despite this com- 
petitive behavior, the FE instability is found to coexist in 
the zero-temperature AFD phase. Once again, quantum 
zero-point fluctuations must be invoked to explain the 
experimental observation that the AFD phase is stable 
against FE distortions. 
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